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Abstract  -  A  theory  is  presented  on  radiation  effects  in  “cellular”  devices 
comprised  of  small  particles  such  as  nanocrystals  or  quantum  dots.  The 
theory  explains  the  surprising  discovery  that  the  photoluminscence  of 
quantum  dot  devices  can  be  significantly  more  radiation-tolerant  than 
bulk  or  quantum  well-based  photodiodes. 
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I.  INTRODUCTION 

In  the  field  of  radiation  effects,  reduced  dimension  has  been  associated  with  increased 
sensitivity  to  radiation  damage.  High  temperature  superconductors  (HTSs),  for  example,  are 
several  orders  of  magnitude  more  sensitive  to  displacement  damage  than  their  low  temperature 
counterparts  because  Cooper  pairs  in  HTSs  are  constrained  to  lie  in  two-dimensional  (2-D)  Cu-0 
planes,  and  when  paired  carriers  scatter  from  radiation-induced  defects  they  can  be  ejected  from 
the  planes  and  hence  from  the  superconducting  state  with  high  efficiency.1  In  contrast,  scattering 
from  radiation-induced  defects  in  isotropic  low  temperature  superconductors  does  not  lead  to  the 
decoupling  of  electron  pairs  because  all  directions  are  essentially  equivalent.  Similarly,  in  high 
electron  mobility  transistors  (HEMTs)  carriers  travel  from  source  to  drain  via  a  2-D  electron  gas 
(2DEG).  Although  radiation-induced  defects  outside  the  2DEG  do  affect  the  operating 
parameters  of  HEMTs,  the  main  displacement  damage  effect  occurs  in  the  2DEG  region,  where 
defects  can  easily  scatter  carriers  out  of  the  transport  state/  Given  these  examples,  it  is 
surprising  to  learn  that  in  irradiated  bulk  (3-D),  quantum  well  (2-D)  and  quasi-zero  dimensional 
quantum  dot  photodiodes  radiation  tolerance  increases  as  dimensionality  is  reduced.  If  this 
effect  could  be  exploited,  the  radiation  tolerance  of  numerous  device  types  susceptible  to 
radiation  damage  could  be  improved.  Unfortunately,  although  a  qualitative  explanation  for  this 
phenomenon  exists,  there  is  as  yet  no  quantitative  theory  to  describe  how  some  low-dimensional 
devices  have  high  radiation  tolerance.  Following  a  brief  review  of  the  available  experimental 
results,  such  a  theory  is  presented  here. 

II.  EXPERIMENTAL  RESULTS 

Results  of  displacement  damage  experiments  have  been  reported  for  a  number  of  bulk, 
quantum  well  (QW)  and  quantum  dot  (QD)  devices,  including  InxGai_xAs/  GaAs  and  GaAs/InAs 
QWs,4"12  and  GaAs/InAs,  InxGai_xAs/GaAs  and  TnAl As/AlGaAs  QDs.4'11 


Generally,  the 
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optoelectronic  performance  of  devices  decreases  with  increasing  particle  fluence  O.  For 
example,  GaAs  pn  junction  bulk  diodes  are  about  30  times  more  sensitive  to  displacement 
damage  than  GaAs/AlGaAs  QW  diodes;  “  and  QW  devices  are  10-100  times  more  sensitive  than 
QD  devices.5’7’8'12 

The  effect  of  radiation  damage  on  the  photoluminescence  P  of  bulk,  QW  and  QD  devices 
is  illustrated  in  Table  1  for  various  materials  systems,  incident  particles  and  particle  energies. 
Normalized  values  of  the  change  in  P,  A P/Pa  are  shown  in  column  7.  To  obtain  a  basis  for 
comparing  different  values  of  A P/P0  the  actual  particle  fluences  (column  5)  were  converted  to 
equivalent  fluences  Oeq  of  2  MeV  protons  (column  6)  by  using  a  standard  conversion  relation 
involving  the  nonionizing  energy  loss  (NIEL,  the  density-normalized  rate  at  which  the  incident 
particle  energy  is  converted  to  atomic  displacements  ).  Although  the  equivalent  fluences  listed 
in  Table  1  are  not  identical,  it  is  clear  that  in  terms  of  photo  luminescence  QD  devices  are 
significantly  more  disorder-tolerant  than  QW  or  bulk  devices.  Similar  trends  have  been  observed 
for  absorption  coefficients  and  electroluminescence  in  QDs  and  QWs. 

The  currently-accepted  explanation  for  increasing  radiation  tolerance  with  decreasing 
dimensionality  is  that  carriers  in  bulk  devices  can  move  in  three  dimensions,  and  so  are 
susceptible  to  scattering,  trapping  and  recombination  at  all  radiation-induced  defects  within  a 
given  volume.  Carriers  in  QWs  are  confined  to  two  dimensions  and  so  are  susceptible  only  to 
those  defects  lying  on  a  disk,  while  carriers  in  QDs  are  dimensionally  constrained  and  so  are  not 
susceptible  to  radiation-induced  defects  outside  their  own  locales.  Within  this  scenario,  carriers 
in  QDs  remain  unaffected  by  radiation  damage  in  other  “cells”. 

The  theory  presented  below  on  radiation  damage  effects  in  charge-localized  systems 
bears  a  resemblance  to  Target  theory,  which  was  developed  in  the  1940s  to  describe  survival 
curves  obtained  by  irradiating  living  tissue  with  ionizing  radiation.14  However,  the  present  Cell 
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theory  allows  for  a  distribution  of  cell  sizes  and  responses,  while  in  standard  Target  theory  all 
unirradiated  cells  are  assumed  to  be  the  same.  Also,  Cell  theory  addresses  displacement  damage 
while  Target  theory  was  developed  to  describe  ionization  damage. 

III.  THEORY 

We  define  a  cell  as  one  of  many  initially  similar  (but  not  necessarily  identical) 
components  which  contribute  piecewise  to  the  function  of  a  device,  so  that  under  the  application 
of  an  external  stimulus  cell  i  exhibits  a  response  r,  in  such  a  way  that  the  device  a  response  R  is 

R  =  Yjri  ’  (1) 

i=i 


where  the  number  of  cells,  D,  is  expected  to  be  large.  For  instance,  a  QD  cellular  device  having 
an  area  A  =  50x50  pm"  and  an  areal  dot  density  pc  =  10  /cm"  (a  commonly  encountered  value) 
would  contain  2.5x1 05  dots.  For  devices  having  more  than  one  layer  of  cells  it  is  convenient  to 
rewrite  Eq.  (1)  as 


(2) 


where  L  is  the  number  of  layers,  D/  is  the  number  of  cells  in  layer  /  and  r$  is  the  response  of  cell 
i  in  layer  /. 

The  problem  at  hand  is  to  determine  how  Eq.  (2)  evolves  as  particle  irradiation  creates 
displacement  damage  in  a  device.  Imposing  a  few  initial  conditions  is  helpful  in  this  regard  (but 
not  necessary;  the  conditions  can  be  relaxed  in  a  more  general  fonnalism).  First,  incident 
particles  are  required  to  have  sufficient  energy  that  they  pass  through  all  cells  without  significant 
energy  loss,  thus  creating  a  uniform  damage  profile.  Second,  incident  angles  are  not  allowed  to 
be  so  steep  that  an  incident  particle  is  likely  to  strike  two  cells  in  the  same  layer.  Third,  layers  in 
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multilayer  devices  must  be  sufficiently  similar  in  composition  that  an  imaginary  average  layer 
can  be  defined  in  such  a  way  that  the  individual  real  layers  are  distributed  smoothly  if  not 
symmetrically  about  a  mean.  Then,  the  term  m  in  Eq.  (2)  can  be  rewritten  as 


rn  =  ria  +  Am  ,  (3) 

where  ria  is  the  the  value  of  r;  in  the  average  layer  and  A m  is  the  deviation  of  the  real  ru  from  that 
average.  Substituting  Eq.  (3)  into  Eq.  (2)  gives 

I  O' \a  +  A r„ )  =  T  Z  ru,  +  X  Z  Arn  '  (4) 

1=1  i= 1  ;=1  1=1  i= 1 

Since  the  individual  layer  structures  are  required  to  be  smoothly  distributed  about  a  mean,  the 
double  summation  over  the  deviations  Am  tends  to  zero  because  cumulative  deviations  tend  to 
cancel  for  large  D.  Hence, 


*  =  •  <5) 

i= l 

The  collapse  of  the  summation  over  /  indicates  that  the  effect  of  irradiation  need  only  be  studied 
for  a  single  layer,  albeit  an  imaginary  one,  then  applied  to  all  layers.  For  convenience,  then,  ria  is 
subsequently  shortened  to  rt,  and  it  is  assumed  without  loss  of  generality  that  L  =  1 . 

Because  by  definition  all  cells  must  be  initially  similar,  their  initial  responses  r,  should 
map  a  smooth  or  at  least  piecewise  continuous  distribution  about  a  mean  value  rc  =  RID. 
Similarly,  the  cross  sectional  areas  A,  of  cells  as  seen  by  an  incident  particle  should  be  distributed 
about  a  mean  value  Ac.  Then,  if  the  device  area  is  A,  the  areal  density  of  cells  pc  is 
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Pc  =  D/A  ,  (6) 

and  the  coverage  C,  of  cells  is 

5  'tIaAIO.  +  M)  ,  (7) 

71  |=1  71  1=1 

where  Ad,  is  the  deviation  of  d,  from  dc.  But  due  to  the  distribution  about  Ac  the  summation 
over  Ad,  tends  to  zero  as  well,  so  Eq.  (7)  reduces  to 


C,  =  D  Ac  /A  =  Ac  pc  .  (8) 

If  the  2.5x10s  QDs  of  the  above  example  have  diameters  near  25  nm  (again,  a  common  value), 
they  cover  about  123  pm  of  the  2500  pm  device  and  result  in  a  coverage  C  ~  0.05. 

Now,  if  an  energetic  particle  impinges  on  a  device,  the  probability  that  it  strikes  any  cell 
is  C,  and  the  probability  that  it  strikes  a  specific  cell  is  C/D.  Thus  for  a  particle  fluence  <E>  the  total 
number  of  incident  particles  is  A  =  Od  and  the  most  likely  number  of  cells  hit  is  n  =  QAC,.  As  an 
initial  exercise  it  is  assumed  that  once  a  cell  is  hit  it  ceases  to  function  (i.e.,  r,  — >•  0).  It  is  also 
initially  assumed  that  ®d  C  «  D,  which  means  that  the  probability  of  a  cell  being  hit  more  than 
once  is  negligible.  Then,  if  the  initial  device  response  is  R0,  the  response  Rn  after  n  =  <E>d<^  cells 
have  been  hit  is 


R  =R  — 'V(r„  -Ar  )  =  R  -nr 

n  o  v  c  j  J  o  c 

7  =  1 


Z  Ar 


7=1 


(9) 
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The  term  A rj  represents  the  initial  deviation  of  the  response  of  cell  j  from  the  average.  In 
experiments,  values  of  N  are  determined  by  the  ion  beam  current,  the  beam  area  and  the  exposure 
time,  and  are  generally  large.  For  the  QDs  and  QWs  of  Table  1  the  maximum  equivalent 
fluences  range  from  about  10  -  10  particles/cm“.  Assuming  a  minimum  fluence  of  10'  times 
the  maximum,  the  previously-mentioned  2500  pm  device  would  receive  a  minimum  of  between 
105  and  106  hits,  resulting  in  a  minimum  of  about  n  =  5000-50000  cells  hit.  Because  n  is  large 
the  summation  over  A/7  in  Eq.  (9)  also  tends  to  zero.  Thus,  substituting  R0/D  for  rc  and  <PA C  for 
n  yields 


m)=R0 


1- 


D 


(10) 


If  the  second  assumption  above  becomes  invalid  and  9) AC,  is  not  much  less  than  D, 
allowance  must  be  made  for  the  possibility  that  a  cell  is  hit  more  than  once.  Then,  it  is  useful  to 
note  that  the  probability  /;(())  that  a  given  cell  remains  undamaged  after  n  cells  have  been  hit  is 


P(0)  =  (1  -  VDf 


(11) 


In  this  equation,  1  ID  is  the  probability  that  a  specific  cell  is  hit  by  one  of  the  n  incident  particles 
to  strike  a  cell. 

After  N  particles  have  hit  the  device  and  n  of  those  have  hit  cells,  the  remaining  number 
of  undamaged  cells  is  Dp( 0)  and  the  overall  device  response  is  rcDp( 0)  (assuming  again  that  the 
deviations  A 77  cancel  upon  summation).  Substituting  for  rc  and  //(())  then  yields 
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R„=R(0)  =  R0(l  -  l/D)®^  .  (12) 

Expanding  Eq.  (12)  to  lowest  order  in  1  ID  (which  for  the  above  example  is  4x1  O'6)  yields  Eq. 
(10),  as  expected  when  OAC  is  small.  Also,  because  D  is  large,  Eq.  (12)  can  be  approximated  by 

R(O)  =  R0  zMQD  (13) 

for  all  but  the  highest  fluences.  Eqs.  (12)  and  (13)  complete  the  initial  exercise. 

An  equation  analogous  to  Eq.  (13)  is  obtained  when  the  so-called  “single-hit”  model  is 
used  in  Target  theory  to  describe  the  survival  curves  of  living  tissue  exposed  to  ionizing 
radiation.  In  Target  theory  the  next  step  in  refining  the  equivalent  of  Eq.  (13)  is  to  switch  models 
from  the  single-hit  to  the  more  elaborate  “multi-target-single-hit”  or  “single-target-multi-hit” 
models.  The  multi-target  model  is  used  to  distinguish  between  various  components  of  a  cell 
(e.g.,  cell  wall,  nucleus)  while  the  multi-hit  model  is  used  to  describe  cells  which  die  after 
undergoing  a  fixed  number  of  ionization  events.  These  models  have  been  used  to  describe 
ionization  effects  in  both  tissue  and  electronic  memory  cells,15  but  are  inappropriate  for 
describing  displacement  damage  effects  in  QDs  and  nanocrystals  because  a  method  is  required 
for  describing  the  gradual  degradation  of  cells.  Thus,  rather  than  using  the  Target  theory  models, 
we  introduce  an  effectiveness  parameter  E  to  describe  the  particle-induced  damage  done  to  a 
cell.  Without  loss  of  generality,  E  =  0  can  be  chosen  to  represent  an  impact  that  does  not  alter  a 
cell’s  response  and  E  =  1  can  be  chosen  to  represent  cell  destruction. 

The  value  of  E  for  displacement  damage  is  expected  to  vary  among  cells  for  several 
reasons.  First,  an  incident  particle  can  impact  a  cell  at  a  variety  of  locations  and  thereby  cause 
different  amounts  of  damage.  Second,  the  stochastic  nature  of  ion  impacts  results  in  various 


9 


amounts  of  damage  even  for  impacts  at  otherwise  identical  sites.  As  a  result,  values  of  E  will 
range  from  near  zero  for  low  energy,  low  mass  particles  striking  a  cell’s  perimeter  to  a  certain 
particle-  and  material-dependent  maximum,  Emax.  (Values  of  Emax  will  probably  depend  on  the 
incident  particle’s  mass  and  energy  via  the  nonionizing  energy  loss.)  As  mentioned,  even  the 
smallest  typical  fluences  result  in  large  numbers  of  cells  being  hit,  which  means  that  infonnation 
about  the  distribution  of  the  values  of  E  is  likely  to  be  difficult  to  extract  from  experimental  data. 
Thus,  only  an  overall  average  effectiveness  parameter  is  expected  to  be  evident.  For 
convenience,  this  average  value  is  called  E. 

Within  the  above  construct  the  function  of  a  cell  after  being  hit  once  is  r,{  1  -  E). 
Therefore,  after  being  hit  m  times,  a  cell’s  function  is  on  average 

rm,i  =  ri(l-E)m  .  (14) 


If  the  particle  fluence  is  sufficiently  small  that  cells  are  unlikely  to  be  hit  more  than  once, 
then  the  most  likely  number  of  cells  not  hit  is  D  -  OAC.  It  is  straightforward  to  use  Eq.  (14)  and 
a  summation  process  similar  to  that  of  Eq.  (9)  to  show  that 


R(®)  =  K 


l- 


O  AE<Z 

D 


(15) 


It  can  also  be  shown  that  if  there  is  a  likelihood  of  multiple  hits  on  a  given  cell,  Eq.  (15)  becomes 


R(0)  =  R0(  1  -  1  /DfAE: 


(16) 


through  a  sequence  similar  to  the  replacement  of  Eq.  (10)  by  Eqs.  (12)  and  (13). 
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Having  now  defined  the  necessary  parameters  and  derived  equations  applicable  to  several 
special  cases,  a  more  comprehensive  development  is  undertaken.  After  n  =  04 C  cells  have  been 
hit,  the  device  funtion  can  always  be  written  as 

n 

Rn=yZN,<ri>  ,  (17) 

1=0 

where  /V,  is  the  number  of  cells  hit  i  times  and  <r,>  is  the  average  function  of  cells  hit  i  times. 
The  probability  that  a  given  cell  has  not  been  hit  is  given  by  Eq.  (11).  The  probability  p(j)  that  a 
cell  has  been  hit  j  and  only  j  times  can  be  shown  to  be 


p(j)  = 


1  y  n\ 

D)  (/;-./)!./! 


(18) 


The  first  and  second  terms  on  the  right  hand  side  of  Eq.  (18)  represent  the  probability  that  a  cell 
is  not  hit  (n-j)  of  the  n  times  but  is  hit  j  times.  The  third  term  represents  the  number  of  ways  the 
hit/no-hit  sequence  can  occur.  Summing  p(j)  over  all  j  from  0  to  n  shows  that  it  is  normalized  to 
unity.  The  number  of  cells  hit  j  and  only  j  times  is  then  Dp(j). 

Rather  than  now  applying  the  mean  value  theorem  as  was  done  following  Eqs.  (7)  and 
(9),  a  related  statistical  principle  is  invoked:  If  an  arbitrary  subset  is  selected  at  random  from  a 
global  distribution,  the  best  initial  estimate  of  the  subset’s  average  values,  standard  deviations, 
and  so  on  is  that  these  values  are  the  same  as  those  of  the  global  distribution.  For  example,  the 
average  initial  function  of  any  randomly  chosen  subset  of  cells  should  be  rc  =  R0/D.  The  actual 
function  of  a  specific  subset  depends  on  how  many  times  the  cells  in  it  have  been  hit.  Using  Eq. 
(14)  to  describe  a  group  of  cells  hit  j  times  and  Eq.  (18)  for  p(j)  allows  Eq.  (17)  to  be  written  as 
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R 


n 


K  Id- ey 


j= 0 


fl-11 

n~J 

r  n 

Dj 

kDj 

n\ 

(■ n-j)ljl 


(19) 


Extracting  a  factor  of  (1  -  1/D)"  from  the  summation  and  introducing  the  variable  y  =  (\-E)/(D-\) 
allows  Eq.  (19)  to  be  rewritten  as 


R 


n 


=  R 


yJn\ 

(n-j)\j\ 


(20) 


Since  E  must  be  between  zero  and  one,  and  since  D  is  large,  the  value  of  y  will  be  much  less  than 
one.  Hence,  yJ  vanishes  as  j  — >  n  and  Eq.  (20)  is  little  changed  if  the  summation  is  perfonned 
from  0  to  oo  rather  than  from  0  to  n.  In  that  case,  the  value  of  the  summation  is  simply  (1  +  yf , 
so 


\  (  1  ) 

n 

E 

1 - 

1  (i  +  y) 

=  R0 

1 - 

A  Dj 

1 

D 

(21) 


Eq.  (21)  reduces  to  the  trivial  case  for  E  =  0  and  to  Eq.  (12)  for  E  =  1.  In  most  practical  cases 
Eq.  (21)  can  be  approximated  by  Eqs.  (15)  and  (16).  Using  Eq.  (6)  for  pc  and  Eq.  (8)  for  C,  yields 
the  final  equation  for  the  displacement  damage  response  of  cellular  devices: 


R. 


R 


i-® 

D 


Re' 


(22) 


IV.  CONCLUSION 

The  response  of  bulk  optoelectronic  devices  to  displacement  damage  commonly  depends  on 
parameters  such  as  carrier  diffusion  lengths,  recombination  constants,  donor  and  acceptor 
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concentrations,  defect  trap  levels  etc.  For  cellular  devices,  however,  the  predicted  response 
depends  almost  exclusively  on  the  particle  fluence,  the  cell  area  and  the  damage  effectiveness 
parameter.  Thus  cellular  devices  do  not  respond  to  radiation  damage  in  the  same  way  that  bulk 
and  (probably  to  a  lesser  extent)  QW  devices  do. 

The  Cell  theory  presented  here  for  describing  displacement  damage  effects  in  quasi-zero 
dimensional  devices  such  as  QD  photodiodes  represents  the  first  time  that  charge  localization  has 
been  built  into  a  model  of  radiation  effects.  This  theory  can  be  used  to  predict  that  the  radiation 
tolerance  of  cellular  devices  can  be  improved  simply  by  employing  smaller  cells.  Finally,  if  Eq. 
(22)  is  validated  by  experiments,  Cell  theory  could  have  a  significant  impact  on  the  way  future 
radiation-hard  devices  are  designed 

The  author  wishes  to  thank  E.A.  Burke  and  G.P.  Summers  for  helpful  comments.  This  work  was 
sponsored  in  part  by  the  Office  of  Naval  Research. 
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Material 

Device 

Particle 

NIEL 

<t> 

A PIP0 

Ref. 

GaAs  pn 

Bulk 

50  MeV  H+ 

3.45E-3 

5E10 

5.56E9 

0.87 

12 

GaAs/InGaAs 

QW 

1.5  MeV  H+ 

5E-2 

3.4E9 

5.66E9 

~0  * 

10 

GaAs/InGaAs 

QD 

1.5  MeV  H+ 

5E-2 

3.4E9 

5.66E9 

~0  * 

10 

GaAs/lnGaAs 

QW 

2.4  MeV  H+ 

2.5E-2 

1E13 

8.1E12 

0.9925 

8,11 

GaAs/InGaAs 

QD 

2.4  MeV  H+ 

2.5E-2 

2E13 

1.6E13 

0.71 

8,11 

GaAs/InGaAs 

QW 

1.5  MeV  H+ 

5E-2 

6E13 

9.7E13 

0.997 

9,10 

GaAs/InGaAs 

QD 

1.5  MeV  H+ 

5E-2 

6E13 

9.7E13 

0.4-0. 8 

9,10 

GaAs/InGaAs 

QW 

1.5  MeV  H+ 

5E-2 

4E13 

6.4E13 

0.995 

5 

GaAs/InGaAs 

QD 

1.5  MeV  H+ 

5E-2 

4E13 

6.4E13 

0.75 

5 

Table  1.  Radiation-induced  changes  in  the  normalized  photo  luminescence,  A PtP0, 
for  selected  bulk,  QW  and  QD  devices.  Interpolated  data  are  marked  by  an  asterisk 
(*).  Units  of  NIEL  are  MeV-cm  /g.  Units  of  fluence  are  particles/cm".  <t>eq  is  the 
equivalent  fluence  of  2  MeV  protons  as  described  in  the  text.  Despite  variations  in 
<E>eq,  values  of  AP/P0  are  consistently  smallest  for  QDs,  indicating  they  are  the  most 
radiation-tolerant. 
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